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Abstract 

The problem of enumerating dimers on an x A/" net embedded on non-orientable surfaces 
is considered. We solve both the Mobius strip and Klein bottle problems for all M. and N with 
the aid of imaginary dimer weights. The use of imaginary weights simplifies the analysis, and 
as a result we obtain new compact solutions in the form of double products. The compact 
expressions also permit us to establish a general reciprocity theorem. 
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1 Introduction 

A seminal development in modern lattice statistics is the solution of enumerating close-packed 
dimers, or perfect matchings, on a finite Ai x Af simple-quartic net obtained by Kasteleyn Q and 
by Temperley and Fisher more than 40 years ago. In their solutions the simple-quartic net 

is assumed to possess free or periodic boundary conditions In view of the connection with 
the conformal field theory where the boundary conditions play a crucial role, there has been 
considerable renewed interest to consider lattice models on non-orientable surfaces ^ ^, ^, |9|. 
For close-packed dimers the present authors have obtained the generating function for the 
Mobius strip and the Klein bottle for even M and Af. Independently, Tesler |10| has solved the 



problem of perfect matchings on Mobius strips and deduced solutions in terms of a g-analogue of 
the Fibonacci numbers for all {Ai,M}. It turns out that the explicit expression of the generating 
function depends crucially on whether A4 and N" being even or odd, and the analysis differs 
considerably when either A4 or AT is odd. In this paper we consider the general {A4,J\f} problems 
for both the Mobius strip and the Klein bottle by introducing imaginary dimer weights. The use 
of imaginary weights simplifies the analysis and as a result we obtain compact expressions of the 
solutions without the recourse of Fibonacci numbers. The compact expressions of the solutions 
also permit us to establish a reciprocity theorem on the enumeration of dimers. 

2 Summary of results 

For the convenience of references, we first summarize our main results. Details of derivation will 
be presented in subsequent sections. 

Consider an MxJ\f simple-quartic net consisting MM sites arranged in an array of M rows and 
M columns. The net forms a Mobius strip if there is a twisted boundary condition in the horizontal 
direction as shown in Fig. 1, and a Klein bottle if, in addition to the twisted boundary condition, 
there is also a periodic boundary condition in the vertical direction. Let the dimer weights be Zh 
and Zv, respectively, in the horizontal and vertical directions. We are interested in the close- form 
evaluation of the dimer generating function 
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where nh,ny are, respectively, the number of horizontal and vertical dimers, and the summation is 
taken over all close-packed dimer coverings of the net. 

Our results follows: For both A4 and Af even, we have Q 
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where the superscripts refer to the type of the nonorientable surface under consideration. 
For M. even and M odd, we have 



^A4X(^h,^.) = Re 
Z^ui^h^Zy) = Re 
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and for M. odd and even, we have 
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For and N both odd, the generating function is zero. 
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3 The Mobius strip 

It is well-known that there is a one-one correspondence between dimer coverings and terms in 
a Pfaffian defined by the dimer weights. However, since terms in the Pfaffian generally do not 
possess the same sign, the evaluation of the Pfaffian does not necessarily produce the desired dimer 
generating function. The crux of matter is to attach signs, or more generally factors, to the dimer 
weights so that all terms in the Pfaffian have the same sign, and the task is reduced to that of 
evaluating a Pfaffian. 

These tasks were first achieved by Kasteleyn [|| for the simple-quartic dimer lattice with free 
and periodic boundary conditions, who showed how to attach signs to dimer weights and how to 
evaluate the Pfaffian. Soon after the publication of T. T. Wu pointed out that the structure 
of the Pfaffian, and hence its evaluation, can be simplified if a factor i is associated to dimer weights 
in one spatial direction. Indeed, the Wu prescription requires only uniformly directed lattice edges 
with the association of a factor i to dimer weights in the direction in which the number of lattice 
sites is odd. (If the number of lattice sites is even in both directions, then the factor i can be 
associated to dimers in either direction.) For {A^,AA} = {even, odd} for example, one replaces Zh 
by izh- 
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To see that the Wu procedure is correct, one considers a standard dimer covering Cq in which 
the lattice is covered only by parallel dimers with real weights. Then, the two terms in the Pfaffian 
corresponding to Cq and any other dimer covering Ci have the same sign, since the superposition 
polygon produced by Cq and Ci contains an even number of arrows pointing in one direction as 
well as a factor = — where n is a nonnegative integer. Namely, the superposition polygon 



is "clockwise-odd" [11|. This use of imaginary dimer weights is the starting point of our analysis. 
3.1 Ai = even Mobius strip 

For M even and M = even or odd, we write for definiteness M = 2M, M = N where M, N are 
positive integers. We orient lattice edges as shown in Fig. 1(a). For the time being consider more 
generally that the horizontal dimers connecting the first and AA-th columns have weights z and the 
associated generating function 



2M 



m=0 



mrp 



(8) 



where Tm = Tra{zh, Zy) is a multinomial in z^ and z^ with strictly positive coefficients. The desired 
result is then obtained by setting z = Zh- Note that Tq is precisely the generating function with 
free boundary conditions. 
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Figure 1: Edge orientations of Mobius strips with twisted boundary conditions in the horizontal 
direction. A, B, C, D, E denote repeated sites, (a). M = A, Af = 5. (b). M = 5, Af = 4. 



Attach a factor i to all (horizontal) dimer weights z^ . This leads us to consider the antisymmetric 
matrix 

A{z) = izh{FN - F^) I2M + zJn ® {F2M - F^m) + z{Kn + kI) ® J2M (9) 

where In is the N x N identity matrix, F^ is the transpose of F, and F2M, Kn and J2M are 
matrices of the order given by the subscripts. 
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(10) 



Now, the Pfaffian of A{z) gives the correct generating function Tq in the case of z = |11]. For 
general z we have the following result: 
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Theorem: The dimer generating function for the simple- quartic net with a twisted boundary 
condition in the horizontal direction is 



1 

2 

Re 



(1 - i)FiAiizh) + (1 + i)FiAi-izh] 
(1 - i)FiA{izh] 



(11) 



Proof. It is clear that the term in corresponding to the configuration Cq {m = 0) has 
the correct sign. For any other dimer configuration Ci, the superposition of Cq and Ci forms 
superposition polygons containing z edges. We have the following facts which can be readily 
verified: 

(i) The sign of a superposition polygon remains unchanged under deformation of its border 
which leaves Uz, the number of z edges it contains, invariant. 

(ii) Deformations of the border of a superposition polygon can change only by multiples of 
2, and the sign of the superposition polygon reverses whenever changes by 2. 

(iii) Superposition polygons having or 1 z edges have the sign +. 
As a result, we obtain 

PfA(z) ^ ,J\A{z)\ 

= Xq + zXi- z^X2- z'^X^ (12) 

where | • | denotes the determinant of • and 

X„ = r„ + z^To,+A + z^T^+8 + ■■■, a = 0, 1, 2, 3. (13) 

The theorem is now a consequence of the fact Zjl^^{zh, z^; Zh) = Xq + ZhXi + z'^X2 + z\X^. 

Remarks: The theorem holds also for the Klein bottle which, in addition to a twisted boundary 
condition in the horizontal direction, has a periodic boundary condition in the vertical direction 
(see below). For the Mobius strip we have Xi = X^ = when = even. 

It now remains to evaluate PfA(ibz2:/j). To evaluate PfA(ibi2;^) = yJ\A{^iZh)\ we make use 
of the fact that, since -F2M — ^2M commutes with J2M7 the 2MN x 2MN matrix A[z) can be 
diagonalized in the 2M-subspace Q. Introducing the 2M x 2M matrix U whose elements are 

(TT-i\ I ■\m' / 2 . / mm'-JT \ 

we find 

{U~'^{F2M - F2M)U)m,m' = C^i COS (t>m) ^m,m' 

{U~\J2MU)m,m' = * ("1)*'+™ W'> m, m' = 1, 2, 2M, (15) 

where 4)m = rm:/{2M + 1). Thus, we can replace the 2M x 2M matrices in (^ by their respective 
eigenvalues, and express |j4(z)| as a product of the replaced determinants, namely, 

2M 

\A{z)\=^'^^^X{\A^;^\z)\ (16) 

m=l 



m,m' = 1,2,..., 2M, (14) 
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where we have taken out a common factor i from each element of the N x N matrix 
A^;;\z) = 2Z, cos 4>mlN + Zh{FN " F^) + {-1)^'+"'z{Kn + Kl). 



(17) 



The matrix A^^f^ {z) can be evaluated for general z in terms of a ^-analogue of Fibonacci numbers, 
but for our purposes when z = Uzh, the matrix can be diagonalized directly. 

Define the N x N matrix 
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we can rewrite (p^ ) when z = izh as 

A^^\iZh) = 2z^ cos (l)mlN + Zh{TN - tIj, 



(19) 



(-1) 



M+m+l 



{An - 1)-k/2N. 



Now T/vT and commute so they can be diagonalized simultaneously with respective eigenvalues 
e*^" and e~*^", where Q 

Thus, we obtain 
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and result 
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where we have made use of the fact that cos02M+i~m = — cos(/)m, (— 1)^*^+^ 
j2MN _ ^_]^-jA/A^_ thus obtain after taking the square root of (|2^ 



M Af 



pfyi(iz,) = n n 
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-(-1)™, and 
(23) 



The substitution of @) into (^) now yields (|). For N = even the Pfaffian (H) is real and (|) 
reduces to (|2D. There is no such simplification for = odd. 



3.2 M = odd Mobius strip 

For M. = odd and M = even or odd, we write, for definiteness, Ai = 2M + 1, M = N. Since the 
number of rows M is odd, we now attach a factor i to dimers in the vertical direction. Again, 
we assign weights z to horizontal dimers connecting the first and A^-th columns and consider the 
generating function 

2M+1 

Zf,%{zh,z,;z)= z^Tm (24) 

m=0 

^ A similar expression of On given in (17) of Ref. Q contains a typo where M + m + l in the exponent should read 
M + m. This does not alter the results of Ref. H. 
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defined similar to (P). It is readily verified that, with lattice edge orientations shown in Fig. 1(b), all 
terms in Tq have the same sign. It follows that we can use the theorem of the preceding subsection 
where Pf(^) is the Pfaffian of the antisymmetric matrix 



A{z) = Zh{FN - F^) ® hhi+i - iZviN ® (-^2^^+! - F^M+i) + ^ H2M+1, 



with Gn = Kj\[ — KJj and 



2M+1 



Vi 



(25) 
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Apply to (p^) the unitary transformation (|T4|) (with 2M replaced by 2M+ 1). The transformation 
diagonalizes (F2M+1 — F2M+1) in (^) and, in addition, produces 



U ^H2M+lU — (-1)*'^-^2M+1- 

{In®U-^)A{z){In®U) 
Zh{FN - FJj) + 2z^ COS 



Thus, we obtain 
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Writing 
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then the matrix A{z) assumes a form shown below in the case of 2M + 1 = 5, 

zGn\ 
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^(5) , 



Interchanging rows and columns, this matrix can be changed into a block-diagonal form having the 
same determinant, 
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For general M we define the 2N x 2N matrix 
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and use the result 



\B^^'^'HzG^\ = \[l + {-ir]zfr\^, + zf 



(32) 
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to obtain 



1, 



M 

2 TT I ylM/ 



\A{z)\ = \A{z)\ = ^-[l + + zY-X{ \A,^ 

m=l 



(34) 



It therefore remains to evaluate |^2')v'*(-^)l- 

The matrix {z) can be diagonahzed for z = iizh- To proceed, it is convenient to multiply 
from the right by a 2N x 2N matrix (whose determinant is (—1)^) to obtain 



A^;;\iz,) ^ A^^\iz,) 



In <S) 
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where 
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Now Q2N commutes with so they can be diagonalized simultaneously with the respective 
eigenvalues e*^" and e"*^" , where = (2n — 1) /2N, n = 1, 2, • • • 2N. Thus we obtain 



Fn i(-l)^'^+™K^\ 
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and, taking the square root of ( |34D with z replaced by izh 



N/2 M 



Y>iA{^zn) = ^[1 + (-!)-](! + ^)z^J' ■ H 11 N ^^^^ ^^^W^ + 

n=l rra=l 



mvr 



2M + 2 



(38) 



Therefore the generating function vanishes identically for = odd. For = even we substitute 
( p8| ) into dllD and make use of the fact that the two sets sin^(2n — l)7r/2A^ and sin^(4n — 1)tt/2N, 
n = 1, 2, • • • , N/2 are identical. This leads to the result (^). 



4 The Klein bottle 

A Klein bottle is constructed by inserting J\f extra vertical edges with weight z^ to the boundary 
of the Mobius strips of Figs. 1(a) and 1(b), so that there is a periodic boundary condition in 
the vertical direction. The extra edges are oriented upward as the other vertical edges. The 
consideration of the Klein bottle parallels to that of the Mobius strip. Again, we need to consider 
the cases of even and odd M. separately. 



4.1 M = even Klein bottle 

For a 2M x A^ Klein bottle with horizontal edges connecting the first and AA-th row having weights 
z, we have as in (^) the generating function 

2M 

Z^'lui^h, z.- z)=Y. Z^Trr,. (39) 

m=0 
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The dimer weights now generate the antisymmetric matrix 

A^^'^iz) = A{z) - zJn {K2M - K^m) 



(40) 



where A{z) is given by (|9|). Following the same analysis as in the case of the Mobius strip, since 
(i) - (iii) still hold we find the desired dimer generating function again given by the theorem (^) 
or, explicitly. 



(41) 



To evaluate ViA^^'^{z), we note that the 2M x 2M matrices {F2M - K2M - F'^m + ^2m) and 
J2M commute, and can be diagonalized simultaneously by the 2M x 2M matrix W whose elements 
are 
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/4M 
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(42) 



We find 



[W-\F2M - K2M - F^M + K^M)W]m,m' 

{W-'j2MW)m,m' 



(2i sin a. 



-r 



M+m 



6. 
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where am = (2m — l)7r/2M. Diagonalizing A{z) in the 2M-subspace, we obtain 

2M 
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■2MN -Q 
m=l 



(44) 
(45) 



where 

A^^\z) = 2z, sina^lN + zh^Fn - F^) + {-if'+'^ziKN + kJ,). 

This expression is the same as (17) for the Mobius strip, except with the substitution of cos 0m by 
sin Thus, following the same analysis, we obtain 



M N 
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The substitution of (H6|) into (^T|) now gives the result (^). For N = even, the Pfaffian ( |4^ ) is real 
and (0) reduces to 



4.2 M = odd Klein bottle 

For a (2M + 1) x Klein bottle, the inserted vertical edges have dimer weights iz^. The consid- 
eration then parallels that of the preceeding sections. Particularly, the desired dimer generating 
function is also given by (|4l|), but now with 

AK1-(z) = A{z) + izJn ® {K2M+1 - K^M+i)- (47) 
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To evaluate P{A^^^{z), one again applies in the 2M + 1 subspace the unitary transformation which 
diagonahzes F2M+1 - K2M+1 - F2M+1 + ^2M+i- Define 



1 



Using the result 
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1 
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where = (2m-l)7r/(2M+l), m = 1,2, - • • ,2M+1, then the matrix ^(z) = {In'^V'^)A{z){In(. 
V) assumes the form in the case of 2M + 1 = 5, 
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Here Bj!J^^ = z^Q^ + 2zt,sinam. Again, interchanging rows and columns, one changes A{z) into 
the block-diagonal form 



/ -DAT 

-ze^^^G 



-zc-'^^Gn 



\ 



N 



B 



(5) 
N 



B 



(2) 
N 



-ze-'"^GN 



-ze^^^GN B 



(4) 
N 



\ 



Explicitly, for general M, the m-th block is 
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We proceed as in ( ^5| ) by multiplying a 2N x 2A^ matrix whose determinant is (—1)^ from the 
right, and obtain 



A2]J (izh) — A^j^'iizh) 



N 



1 
-1 



2zy sin amhN + Zh{FN - fJj) ® 
2zy sin amhN + Zh{Q2N - QIn)' 



1 
-1 



+ izhGn <8) 



e"*" 



(53) 
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where 



Q2N 
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'Kn 
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(54) 



Now Q2N commutes with <32Af' they can be diagonahzed simultaneously with the respective 
eigenvalues e*^" and e~*^" where j3n = (2n — 1)/2A^, n = 1, 2, • • • , 2N . Thus we obtain 



N 



= n 

Using (|53[) and (|33|), we get 



n=l 



2 ■ 2 - l)7r 



sin 



2A^ 



{2m - l)7r 
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9 9 (2n — l)7r 9 9 (2m — l)ir 
Azf^ sin'' - — ^ + 4z,i sin^ ^ ^ 



2A^ 



2M + 1 



■ (56) 



Thus, the generating function ( pD vanishes identically for = odd. For N = even, we replace as 
before sin (2n - l)Tr/N by sin2(4n - l)7r/A^. The substitution of (|5|) into (g^) now leads to the 
result (^. 



5 A reciprocity theorem 

Using the explicit expression of dimer enumerations on a simple-quartic lattice with free boundaries, 
Stanley [12| has shown that the enumeration expression satisfies a certain reciprocity relation, a 
relation rederived recently by Propp |13| from a combinatorial approach. Here, we show that our 
solutions of dimer enumerations lead to an extension of the reciprocity relation to enumerations on 
cylindrical, toroidal, and nonorientable surfaces. 

We first consider solutions d), (||), (|) and (|) for M = even. Writing 



r^^°^>l, AA) = Zf,%{l, 1), T^'-{M,M) = zKin (1, 1) 



and using the identity [|r 

n-l 

n 

fc=0 



2x cos a + 



2kTT 



+ 1 



n 



X 



2n 



2x" cos (no) + 1 



repeatedly, we can rewrite our solutions (for general z^ and Zy) in the form of 

[iM+l]/2] 
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N/2 

^MMI2 

n=l 



+ (-l)^y-^-i 



Vn+Vn 



ZKln / \ 



MM/2 



[(M+l)/2] 



m=l 



^ 3-(-1)-^ ^.>,^/2^'l.A. 



n [y^ + (-2/n)- 

n=l 



M 



(57) 
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(59) 



(60) 
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where 



^\ — cos ■ . , , yn = G\ — sm ^ — , tm = G\ — sm — (61) 



and G{y) = y + y/y^ + T. Thus, the following reciprocity relations are obtained by inspection: 

T^°^{M,Ar) = T^"^{M, -Af) = exM T^"^{-M - 2,N) 
T^^''{M,J\f) = T^^''{M,-J\f) = ej^MT^^''{-M,M), (62) 

where 

eMM = (-1)^, AA = 2 (mod 4) 

= +1, otherwise. (63) 

There are no reciprocity relations for J\f = odd. 

We have carried out similar analyses for dimer enumerations on a simple-quartic net embedded 



on a cylinder and a torus, using the solutions given in [15, |l|], and have discovered universal rules 
of associating reciprocity relations to specific boundary conditions. Generally, there are 3 different 
boundary conditions, or "matchings", that can be imposed between 2 opposite boundaries of an 
Ai X M net. The conditions can be twisted such as those shown in the horizontal direction in Fig. 
1, periodic such as on a torus, or free which means free standing. To establish the convention we 
shall refer to the boundary condition between the first and the M-th. columns as (the boundary 
condition) in the M direction, and similar that between the first and the A^-th rows as in the 



Ai direction. Then, our findings together with those of Ref. [13| lead to the following theorem 
applicable to all cases: 

Reciprocity theorem: Let T{Ai,J\f) be the number of close-packed dimer configurations (perfect 
matchings) on an Ai x M simple-quartic lattice with free, periodic, or twisted boundary conditions 
in either direction. (The case of twisted boundary conditions in both directions is excluded). If the 
twisted boundary condition, if occurring, is in the Ai (Af ) direction, we restrict to Ai (M ) = even. 
Then, we have 

1. T{Ai,M) = ej\fM T(— 2 — Ai,N) if the boundary condition in the Ai direction is free. 

2. T{Ai,M) = ej\fM T{—Ai,M) if the boundary condition in the Ai direction is periodic or 
twisted. 



6 Summary and Discussions 

We have evaluated the dimer generating function (|l|) for an Ai x J\f simple-quartic net embedded 
on a Mobius strip and a Klein bottle for all A4,M. The results are given by (||) - (|^). Our results 
can also be written in terms of the g-analogue of the Fibonacci numbers J-n{Q) defined by 

oo 

E-^n(^)«"- (64) 



l-qs-s^ 

Using the first line of (|59|), for example, and the identity 



q = X — x~^, (65) 
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one can verify that our results (§) and @ for the Mobius strip are the same as those given by 
Tesler |T0]. Details of the proof which also lead to some new product identities involving the 



Fibonacci numbers will be given elsewhere |16|. We have also deduced a reciprocity theorem for 
the enumeration T(Ai,J\f) of dimers on an 7W x lattice under arbitrary including free, periodic, 
and twisted boundary conditions. 

Finally, we point out that the results (^) - can be put in a compact expression valid for all 
cases as 



[{M+l)/2] AT / //I IN 

n nf2.(-i)^-«dn(^!-^+2x„ 



m=l n=l \ 



(66) 



where [x] is the integral part of x, and 



I Z \ TTITT 

= — COS — for the Moebius strip 

\zhj M + l 

= (^] sin ^^""r/^^ for the Klein bottle. (67) 
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